Transverse Momentum Spectra of J/ip and ip' Mesons 
from Quark Gluon Plasma Hadronization in Nuclear Collisions 
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Recent results on transverse mass spectra of J/tp and t/'' mesons in central Pb+Pb collisions 
at 158 A-GeV are considered. It is shown that those results support a hypothesis of statistical 
production of charmonia at hadronization and suggest the early thermal freeze-out of J/tp and tp' 
mesons. Based on this approach the collective transverse velocity of hadronizing quark gluon plasma 
is estimated to be (vqf) ~ 0.2. Predictions for transverse mass spectra of hidden and open charm 
mesons at SPS and RHIC are discussed. 
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The idealized concepts of chemical (hadron multi- 
plicities) and thermal (hadron momentum distributions) 
freeze-outs were introduced to interpret data on hadron 
production in relativistic nucleus-nucleus (A+A) colli- 
sions 1^. The first experimental results on yields and 
transverse mass (my = ^/Pt ~^ "^^) spectra suggested 
the following scenario: for the most abundant hadron 
species (tt, TV, K, A) the chemical freeze-out, which seems 
to coincide with the hadronization of the quark gluon 
plasma (QGP), is followed by the thermal freeze-out oc- 
curring at a rather late stage of the A-l-A reaction. In this 
letter we discuss whether the new data of NA50 on 
transverse mass spectra of J/tp and ip' mesons produced 
in central Pb-|-Pb collisions at 158 A-GeV are consistent 
with the above picture. Our consideration is based on 
a recent hypothesis Q of statistical production of char- 
monia at hadronization. We further suggest that thermal 
freeze-out of charmonia coincides with the hadronization 
transition. The consequences of this assumptions are in 
agreement with existing data on rriT spectra. They allow 
to extract the transverse flow velocity of the hadronizing 
QGP and lead to new predictions which can be tested by 
future measurements. 

Rescattering among partons and, in the late stage of 
the reaction process, hadrons created in relativistic A-l-A 
collisions should cause local thermalization and the de- 
velopment of transverse collective flow of matter. Thus 
the flnal transverse motion of hadrons can be considered 
as a convolution of the transverse flow velocity of the 
freezing-out matter element with the thermal motion of 
the hadrons in the rest frame of this element. The result- 
ing niT spectrum has approximately exponential shape: 



rriT drriT 



(1) 



parameters, respectively. The T* parameter is related 
to the thermal freeze-out temperature Tf and the mean 
transverse flow velocity {v^). In the nonrelativistic ap- 
proximation (tot < 2to) contributions from thermal and 
collective particle motions can be separated and the in- 
verse slope parameter of the observed hadron spectrum 
can be expressed as Q : 



T* = Tf + - ■ TO ■ 

TT 
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Note that factor 2/7r appears in Eq.(||) due to the cylin- 
drical geometry expected in high energy A-l-A collisions 
(see for details). In the case of the simultaneous ther- 
mal freeze-out of all hadrons the inverse slope parameter 
T* should follow the linear dependence on m given by Eq. 
(^). In Fig. 1 we present a compilation of inverse slope 
parameters measured for various hadron species in cen- 
tral Pb-f-Pb collisions at 158 A-GeV. It is observed that 
the simplistic hypothesis of common thermal freeze-out 
of all hadrons is not supported by the data. 

The freeze-out condition for pions can be determined 
independently of the other hadron species from the re- 
sults on two pion correlations. This procedure leads to 
the following values of the pion thermal freeze-out pa- 
rameters [S|: 



Tf 



120 ± 12 MeV 



(4) 



0.5 ±0.12 



(3) 



where C and T* are a normalization and an inverse slope 



In Fig. 1 the upper solid line indicates the dependence 
given by Eq. (^ for the pion freeze-out parameters. The 
values of T* for and p are close to those given by 
this line (^,^. The results on T* for other hadrons are 
in general below the 'pion freeze-out line'. This may be 
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FIG. 1. The inverse slope parameter as a function of 
the particle mass for central Pb+Pb collisions at 158 A-GeV. 
The results for the following hadrons are compiled from: vr 
mesons, K mesons and protons [6] (closed triangles up), p 
and uo mesons [7] (open square), (j) mesons [8] (open trian- 
gles up) and [7] (open square) , A hyperons [9] (closed triangle 
down), H hyperons [10] (open triangles down), Q. hyperons [11] 
(closed diamonds) , J/tj) mesons [2] (closed square) . The filled 
circles are predictions of the model for ^' and D mesons. The 
upper solid line corresponds to Eq. (2) with the freeze-out 
parameters (3) extracted from the two pion correlation data 
[5]. The lower solid line is given by Eq. (5) with parameters 
Th = 175 MeV and {v") = 0.19 which correspond to the 
QGP hadronization. 

interpreted in the following way. For hadrons with "low" 
masses and large interaction cross sections the thermal 
freeze-out takes place at the late stage of the expansion, 
i.e., at large (v:^) and small Tj given by Eq. (3). Heavy 
hadrons (due to large mass and smaller cross section) are 
expected to decouple from the system early thus leading 
to the smaller values of T* than those expected from the 
'pion freeze-out line' ^. This is seen from the mr-slopes 
of the hyperons and most clearly can be illustrated by a 
small T* value of the f2 hyperon shown in Fig. 1. This ex- 
planation is in the line with the typical hydrodynamical 
calculations (see, for instance, [|l2|) - the inverse slope pa- 
rameter T* of the hadron is smaller, if it decouples early. 



i.e., at higher temperatures Tf, but smaller velocity {v^). 
A more careful analysis requires the combination of the 
hydrodynamic approach at the early stage of the expan- 
sion with cascade model calculations at the latest stage 

ilia. 

The yields of produced light hadrons are surprisingly 
well reproduced within the statistical approach to parti- 
cle production. The temperature parameter extracted 
from the fit to the multiplicity data is found to be 
Th = 175±10 MeV [||-ll|]. This chemical freeze-out 
temperature appears to be quite close to the expected 
temperature of the transition between the hadron gas and 
the quark giuon plasma. This fact suggests the possibil- 
ity to ascribe the observed statistical properties of hadron 
yield systematics at high energies to the statistical nature 
of the hadronization process, i.e., the hadrons are "born" 
at hadronization temperature Th into the state of chem- 
ical equilibrium and their multiplicities are frozen-out 
afterwards. 

Within this approach one can make a rough estimate of 
a lower limit for the value of the measured inverse slope 
parameters: 



-MIN 



175 MeV. 



(4) 



This is done under two extreme assumptions: there is 
no collective transverse flow of hadronizing matter and 
there is no rescattering between produced hadrons. The 
measured values of T* parameter for all hadrons satisfy 
the condition T* > T^,jjj^ showing that the approach 
sketched above leads to self-consistent results in the light 
hadron sector. 

A very different approach is traditionally used in 
modeling the production process of heavy hadrons like 
quarkonia. It is based on the assumption that quarkonia 
are created significantly prior the hadronization of QGP. 
This in general implies that the systematics of quarko- 
nium production should be different than the one estab- 
lished for light hadrons. However recently it was found 
that multiplicity of the J/ mesons in nuclear collisions 
at high energies follows dependences well known from 
light hadron data [|l8|,^ and that it is consistent with 
the hypothesis of statistical production of J /ip mesons at 
hadronization In particular the data on J/ip and ip' 
yields in central Pb-|-Pb collisions at 158 A-GeV are con- 
sistent with the predictions of statistical models 20 -2^] 
for a typical values of Th = 175 MeV extracted from light 
hadron systematics. 



^For pions, in contrast to other hadrons, the relativistic ef- 
fects neglected in Eq. (^) are important. This is a main reason 
why the measured pion inverse slope parameter is above the 
thermal pion freeze-out line in Fig. 1. 
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The new hypothesis of statistical J/ijj production 
at hadronization can be further tested using data on 
Tot spectra. As it follows from the above discussion 
one expects for J /ij) mesons an exponential shape of 
the rriT distribution with the inverse slope parameter 
T*{J/'ip) > Th- Recently pubhshed experimental results 
of NA50 1^ on transverse mass spectra of J/ip mesons 
in central Pb+Pb collisions at 158 A-GeV confirm this 
expectation: the shape of the spectrum is approxi- 
mately exponential with the fitted inverse slope param- 
eter T*{J/tJj) = 245 ± 5 MeV. The measured T* {.]/%!)) 
value is significantly smaller than that expected on the 
base of the pion freeze-out line (T* = 610 MeV). 

The 'low' value of T* for J/^ suggests its rather early 
thermal freeze-out. One may argue that this is due to the 
large mass and low interaction cross section of the J /"ip 
meson. We postulate therefore that the thermal freeze- 
out of J /'ip coincides with the hadronization of QGP, i.e., 
that the J /ip meson does not participate in the hadronic 
rescattering after hadronization. It is, however, natu- 
ral to expect that there is a significant collective trans- 
verse fiow of hadronizing QGP developed at the early 
stage of partonic rescattering. Consequently the inverse 
slope parameter of J/ip meson as well as all other hadrons 
for which chemical and thermal freeze-outs coincide with 
hadronization can be expressed as: 



dependence (||) of the apparent inverse slope with ap- 
proximately the same value of Th = 175 MeV, but with 



(5) 



where {v^ ) is the mean transverse fiow velocity of the 
QGP at the hadronization. Assuming Th = 175 MeV 
and using measured value of T*{J/tp) = 245 MeV we 
find from Eq.(||): (v!j!) = 0.19. As expected the ob- 
tained transverse flow velocity of QGP at hadronization 
is significantly smaller than the transverse fiow velocity 
of pions (« 0.5). The linear TO-dependence of T^ (|^) 
is shown in Fig. 1 by the lower solid line. Within the 
approach discussed here, Eq. (||) can be used to obtain a 
next estimate of the lower limit of the measured inverse 
slope parameters for all hadrons. In fact the values of 
the parameter T* for all light hadrons are higher than 
T^{m). The recent results ^ on the tot spectra of the 
-0' meson indicate that T*{ij') « T^it/j') = 258 ±5 MeV, 
which suggests that also the tp' meson (like J/ip) does 
not participate in the hadronic rescattering. 

One may expect that the thermal freeze-out may coin- 
cide with hadronization also for the D meson. Under this 
assumption we calculate the value of the apparent tem- 
perature for the D meson: T*{D) ^ T^{D) ^ 217 MeV. 
Note that this result is significantly lower than the pre- 
dictions of the model of Ref. although in the 
same value of the hadronization temperature was used. 

Another important issue is the production of open and 
hidden charm particles in A-l-A collisions at RHIC ener- 
gies. One can expect stronger transverse collective flow 
effects than at the SPS. This will lead to a linear mass 



a larger value of 



It is necessary to mention that a 



recent analysis | |25| | of the particle number ratios of the 
RHIC data leads to the above value Th = 175 MeV of 
the chemical freeze-out temperature. The preliminary 
RHIC data of hadron tot spectra measured by STAR 
|6| and PHENIX ^ collaborations support the fact of 
a larger value of the transverse velocity. However, the 
measurements were made for tt^, ,p, and p only, and, 
therefore, this corresponds to the late thermal freeze-out 
stage. 

As long as the RHIC data on the inverse slope pa- 
rameters for the open and hidden charm mesons are 
absent it is interesting to compare the different model 
predictions. Assuming Th = 175 MeV we need the 
value of (w-^ ) to estimate the inverse slope parame- 
ters for charmed hadrons in our approach. The hydro- 
dynamic calculations of Ref. predict the value of 
{v^) = 0.30 at the hadronization in Au-fAu collisions 
at RHIC. This leads to an increase of the inverse slopes 
of charmed hadrons at RHIC in comparison to those val- 
ues at SPS, e.g., T*{J/iP) ^ T^{J/^) = 350 MeV and 
T*{D) ^ T;j{D) = 280 MeV. 

Note that a signiflcantly larger value of 
T*{D) ^ 380 MeV in Au-hAu collisions at RHIC would 
be obtained, if the thermal freeze-out of D mesons hap- 
pens at temperature of 130 MeV. On the other hand, 
the calculations done with the PYTHIA generator for 
p + p collisions show an even higher value of the inverse 
slope for D mesons T*{D) ~ 440 MeV (l|. In fact, 
the situation is even more uncertain as the dynamical 
transport calculations of the HSD model (Hadron String 
Dynamics) made for RHIC energy |28| predict a rather 
small inverse slope of about 225 MeV for all particle 
species. Therefore, the T*{m) systematics at RHIC for 
Tot — to < 1 GeV will provide an opportunity to find 
whether the statistical hadronization works, whether a 
thermal freeze-out of charmed particles happens simulta- 
neously with hadronization and whether the transverse 
collective fiow at hadronization of the QGP is stronger 
at RHIC than at SPS. 

In summary, recent results on transverse mass spec- 
tra of J/tp and -0' mesons in central Pb-|-Pb collisions at 
158 A-GeV are considered. It is shown that these data 
support the hypothesis of the statistical production of 
charmonia at hadronization and suggest a simultaneous 
hadronization and the thermal freeze-out for J/ip and tp' 
mesons. Based on this approach the collective transverse 
velocity of hadronizing quark gluon plasma is estimated 
to be (w^) « 0.2. Prediction for transverse mass spectra 
of hidden and open charm mesons at SPS and RHIC are 
discussed. 
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